Quasiparticle transitions in a charge-phase qubit probed by rf-osciUations 
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We investigated transitions in an Al charge-phase qubit of SQUID-configuration which was in- 
ductively coupled to an rf tank circuit that made it possible to read out the state by measuring the 
Josephson inductance of the qubit. Depending on the flux and charge bias and on the amplitude 
of the rf-oscillations, we probed either the ground state or a dynamic change of the qubit states 
which we attributed to stochastic single quasiparticle tunneling onto and off the qubit island, in- 
volving an exchange of energy with the qubit. Within the scope of this model, a selection rule for 
quasiparticle-induced transitions in the qubit is discussed. 

PACS numbers: 74.50.-|-r, 85.25.Cp, 73.40.Gk 



Superconducting circuits based on small Josephson 
junctions and enabling the tunneling of single Cooper 
pairs offer great opportunities for electronics and quan- 
tum computing. The class of the siiigle Cooper pair de- 
vices comprises Bloch transistors P, 0, [3 , quantum elec- 
trometers [1, H, 2e-pumps (tJ, etc. A prerequisite for 
the regular operation of these devices is that no quasipar- 
ticle (QP) transitions occur and that the so-called even- 
parity state of the small superconducting island is thus 
maintained. If this requirement is not fulfilled, even a 
single QP can instantly change the island charge, leading 
to a shifting of the operation point. 

There are several ways of observing the even-odd states 
of the Bloch transistor island experimentally. These in- 
clude: (i) the measuring of the gate dependence of the 
Josephson switching current (see, e.g^Ref.jl]), (ii) the 
monitoring of the supercurrent peak [3, Q , and (iii) the 
measuring of the island charge q in the zero-current- 
biased transistor (i.e. Cooper pair box) by means of a 
capacitively coupled SET electrometer [l^, [ll|, [l^ . 
Recently, Naaman and Aumentado have investigated 
the parity states of a current-biased Bloch transistor in- 
cluded in a resonance LC-circuit driven by a 500 MHz 
harmonic signal. In this circuit, the charge-dependent 
parameter was the Josephson inductance L]{q) of the 
Bloch transistor in the ground state. 

In this paper, we will address the problem of single 
QPs tunneling in Josephson qubit circuits. The opera- 
tion of charge qubits is, by nature, sensitive to the 
incoherent tunneling of unpaired electrons because they 
change the charge state of the system instantaneously 
and stochastically [l^, [l^ . This makes the setting of the 
optimal working point of the qubit difficult and increases 
the decoherence [H, 

Our circuit comprises a Bloch transistor which is in- 
cluded in a superconducting ring inductively coupled to 
an rf-driven tank circuit [T?, IS!. [l9l [20| . In such an rf- 
SQUID configuration, the Bloch transistor is galvanically 
decoupled from the measuring circuit, which in general 
leads to a reduction in the density of non-equilibrium QPs 
that are able to enter the island. The effect of the QP 
transitions manifested itself in excitations of this qubit 
and was detected by measuring the change in the Joseph- 



son inductance value Lj (g) . We will demonstrate that in 
the general survey, the dynamics of our circuit can be 
described in a similar way as in the model of Lutchyn et 
al. [2l| , which has been published recently and addresses 
the problem of the energy relaxation in the charge qubit 
caused by a single QP. In contrast to Ref. [2111, however, 
we will focus in this paper on the reverse process, i.e. on 
the energy transfer from the QP to the qubit, leading to 
an excitation of the qubit. Moreover, we will derive a se- 
lection rule showing that the rates of the QP transitions 
strongly depend on the initial qubit state. 

In our Al qubit circuit (see Fig.[T] (a)), the super- 
conducting loop which closes the Bloch transistor has 
a small geometrical inductance L Lj. Thus, the to- 
tal inductance of the closed loop is determined mostly 
by the transistor's Josephson inductance. The transis- 
tor is operated as a Cooper pair box (qubit) whose dis- 
tinct quantum states - which have the ground state en- 
ergy Eo and the excited state energy Ei (see Fig. 1(b)) 
- are associated with different Bloch-bands of the sys- 
tem 22| . The quantum states of the transistor also de- 
pend on the phase coordinate ipdc determined by the 
external magnetic flux <I>dc that is applied to the loop 
fdc — 27r<I>dc/*I'oj with $0 being the flux quantum. The 
Josephson inductance Lj is related to the local curva- 
ture of the corresponding energy surface En{q,ip), i.e. 



l/Lj{n,q,ip) cx E^{q,ip) / dip^ ,n 



and 1. Thus, 



the qubit eigenstates can be identified by means of the 
Josephson inductance of the transistor which is probed by 
small rf-oscillations in the loop, with the resulting phase 
Lp{t) — ipdc + a sin (2-Kft) and with a being proportional 
to the amplitude of the rf-oscillations in the tank circuit 
induced by an rf-driving current. The black arcs on the 
energy surfaces in Fig. 1 (b) illustrate the rf-oscillations 
of the Josephson phase, which in turn, show the curvature 
of the band profile and are detected by the tank circuit 
coupled to the qubit. The drive frequency / is close to 
the bare resonance frequency /o ~ 77 MHz of our Nb- 
tank circuit which has a quality factor oi Q k, 370. The 
symmetric double loop is coupled - through the mutual 
inductance M - to the coil of the double spiral tank circuit 
formed as a planar gradiometer. The coupling coefficient 
IS k = M/y/LL^ « 0.4, with M « 3.8 nH, L w 0.7 nH 
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FIG. 1: (a) Diagram of the measurement set-up. The core el- 
ement is a double Josephson junction (the two crossed boxes) 
with a capacitive gate coupled to its island, i.e. the Bloch 
transistor, embedded in a macroscopic superconducting loop. 
The loop is inductively coupled to an rf-driven tank circuit 
which is capacitively coupled to a cold preamplifier, (b) En- 
ergy band diagram of the circuit, calculated for the param- 
eters found in the experiment. The thick-line arcs show the 
variations of the phase ip caused by oscillations in the tank cir- 
cuit. The arrows denote the quasiparticle-induced transition 
between the qubit states (the dashed arrows indicate transi- 
tions suppressed due to the destructive interference effect, see 
text below), (c) Phase shift a in the tank circuit, measured as 
a function of the dc- value of the Josephson phase </5dc and the 
island charge q for an amplitude a — 0.56 rad of rf-oscillations 
of the phase. 



and the tank circuit inductance Lt ~ 150 nH (see Ref. 
[20l | , where the behavior of an all-Nb circuit of similar de- 
sign was investigated). Due to the coupUng to the qubit, 
the effective inductance L^s of the resonance circuit is 
equal to Lt — M'^LJ^{n, q, ip). 

The Bloch transistor, the loop and the tank circuit 
were fabricated by electron beam lithography on the 
same chip. The tank circuit inductor was fabricated on 
the basis of Nb technology [13, Ull , the qubit loop and 
the Bloch transistor by means of the two-angle Al shadow 
evaporation technique. No special precautions for the 




suppression of QP-poisoning of the island - such as, for 
example, the engineering of a barrier-like gap profile with 
the island gap value being greater than the electrode gap 
2J,l25|, or the implementation of normal- metal QP 
in the outer electrodes - were taken. The critical 
currents of the single junctions were approx. 25 nA, with 
the corresponding value of 45 fieV for the average Joseph- 
son coupling energy Ejo = {E,n + E,-;2)/2 (« En « E}2, 
whereby the corresponding values of the single junc- 
tions of the transistor yie ld the effective Josephson cou - 
pling energy E]{ip) = y^E]^ + Ej^ + 2EjiEj2 cos(^dc)- 
The charging energy Eq of the transistor island has a 
value of 110 /^eV, so that both energies Ec and i?jo 
are smaller than the value of the superconducting gap 
in Al films, Aai ~ 210 //eV. These values of Ejq and 
Ec were taken from a fitting of the shape of the ground 
state extracted from rf-measurements with finite ampli- 
tude of the Josephson phase oscillations, see Ref.f2^ for 
details. Moreover, the obtained data agreed well with 
the corresponding Ambegaokar-Baratoff and Coulomb- 
blockade values of similar test transistors, extracted from 
dc-measurements of their TV characteristics. 

Unfortunately, due to the technical conditions (poor 
signal-to-noise ratio for a small amplitude of the phase 
oscillations), it was not possible to determine the Joseph- 
son energy difi'erence Ai?,j = \Eji — E}2\ = E}{ipdc — 
tt) <gi Ejq. An analysis of the fitted curves a{q) shows 
that their shape for the smallest amplitude a = 0.28 of 
the phase oscillations does not depend on the asymme- 
try factor AEj / Ejq if this ratio is smaller than 0.2. From 
this we conclude that the asymmetry factor in our sample 
is in the range AEj/Ejo < 0.2. 

In our experiment, which was carried out in a dilution 
refrigerator at a base temperature of 20 mK, we mea- 
sured the phase angle a between the driving signal Irt 
and the rf voltage K-f on the tank circuit. From the 
a-dependence one can deduce the Josephson-inductance 
Lj(g, ip) by applying the simple formula 



tana = k Q 



(1) 



which applies when a drive frequency / is equal to the 
bare resonance frequency /q. The measurement of the 
phase shift, a, for different values of the dc phase and the 
gate charge allows the curvature of the energy surface to 
be mapped. Figure[T] (c) shows a periodical dependence 
of a both on ip^c = (27r/$o)Af/dc and q = CqVg + goffsot, 
where we find a 2e-periodic gate modulation. As had 
been expected for the ground-state shape, a plain gate 
modulation dependence with a peak-to-peak amplitude 
of 0.6 rad is to be seen in Fig.[T] (c) for (p^c — 0. As 
can be expected from the theoretical dependence of the 
ground state energy, the gate modulation curve for ip^c — 
TT is inverted with respect to the (pdc — 0-curve, which 
is expressed by the opposite sign of the energy surface 
curvature. A closer look at the gate-dependence curve 
measured at (p^c = tt reveals, however, a peculiar shape. 
Before reaching the zenith of the curve, the phase angle 
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a starts to rise sharply and remains in a broad range 
around q — e ed a level that is even higher than that 
for g = (this dependence is also shown in Fig. 2(a) by 
symbols). 

Om- theoretical understanding of this behavior bases 
on the assumption of a statistical mixture of the differ- 
ent quantum states of the qubit. Here, at the degeneracy 
point {q = e and tpdc = tt), the states considered are the 
ground state A and the excited state B, both with even 
parity of the island, and the ground state C and the ex- 
cited state D, both with odd parity, as shown by the black 
arcs in Fig. 1 (b). In our opinion, the excited odd-parity 
state D does not play any significant role, as its exci- 
tation energy is too great (i.e. around 3Ec) and hence 
much larger than the energy gap AEj between A and B. 
We can also rule out a notable contribution of the ground 
state C with odd parity, as the admixture of the values 
of a{q) corresponding to this state (indicated by the cor- 
responding black arcs in Fig. 1 (b)) with small negative 
curvature cannot yield the above-mentioned overshooting 
of the experimental data in the vicinity oi q = e. Like- 
wise, we are of the opinion that the contribution of the 
odd state is small due to the presumably short lifetime 
of a QP in the island, see, e.g., Ref.iUlll. Therefore, 
we modelled the experimental data as a mixture of the 
ground state A and the excited state B, with the signs 
opposite to each other. This was done by replacing, in 
Eq. Il]), the Josephson inductance with 



(a) -0.5 



Lj ^{q) Woiq)Lj ^(g)|„=o -t- wi{q)Lj \q)\. 



(2) 



(whereby the weight factors are non-negative occupation 
numbers which obey the relation 'Wo{q) + wi{q) = 1) 
and have been fitted subsequently to the measured a{q)- 
dependence, see Fig. [5] (a). In this way, we were able to re- 
construct the occupation numbers wo,i{q) of the ground 
state and of the excited state, see Fig.[2](b). As a result, 
we found a,t q = ±e an increase in the occupation of the 
upper state up to wi(e) = 0.46, which remained rather 
large in a broad range around this degeneracy point. 

The found steady-state populations wo.i yield the 
ensemble-averaged values of the Josephson inductance 
according to Eq. (2). The ratio of wi to wq reflects the 
competition between the rates of the QP tunneling and 
the energy relaxation of the qubit. Such a relaxation 
also occurs due to the coupling to the environmental de- 
grees of freedom (e.g. flux and gate control lines, ex- 
ternal magnetic field, background charge, etc.). From 
our measurements at (p^c = tt we can conclude that the 
relaxation due to the environmental degrees of freedom 
which reduces wi and increases wq is not dominant, be- 
cause in our case wq ^ wi ~ 0.5. Since we measure with 
our technique only averaged values of the phase shift a, 
we are not able to find absolute values of the relaxation 
rates. Beyond this, an accurate determination of the dif- 
ferent relaxation rates is only possible when the distri- 
bution of the non-equilibrium QPs is known with regard 
to their energy. To explain this mixture effect, we rule 
out a thermal activation of the excited state at the given 
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FIG. 2: (a) Comparison between the experimental (symbols) 
and the calculated (sohd line) gate modulation dependencies 
a{q) for an amplitude of rf-oscillations of a = 0.56 rad and 
a phase bias i^^dc ~ tt. (b) Reconstructed occupation prob- 
abilities of the ground state (dotted line) and of the excited 
state (solid line), yielding the modelled dependence in (a), 
(c) QP-pumping mechanism: in the left panel, an unpaired 
non-equilibrium QP tunnels from the outer electrode onto the 
island of the qubit in the ground state. In the favoured pro- 
cess, the QP tunnels - due to the somewhat larger density of 
states there - back into a lower-energetic state of the outer 
electrode and, hence, transfers energy to the qubit by excit- 
ing it to the upper state (see panel at the lower right). The 
alternative process, i.e. the tunneling of a QP back to the 
outer electrode without the qubit being excited is shown in 
the panel at the upper right. 



base temperature, which is below 100 mK, because at 
such low temperatures such high wi-values can never be 
achieved, i.e. the Boltzmann-factor for the excited state 
is negligible. Besides, the broad range of the gate charge 
where mixture occurs (more than 40% of the gate charge 
period) cannot be achieved by such a weak thermal exci- 
tation. This is shown by the following estimation: for a 
finite detuning from the degeneracy point q — e, the gap 
energy increases roughly like AE — AEc\q/e~ 1| and the 
factor exp (— AE'/fce?^) is therefore strongly suppressed 
at a finite detuning 5q = q/e — 1 = 0.1 for Eq « 110 /xeV 
and T = 100 mK. 

Another possible explanation for the observed mixture 
effect could be Landau-Zener (LZ) tunneling due to the 
rf-drive which leads to a periodic passing of the degen- 
eracy point (at ipAc — and q — e). When we take the 
broad range of the mixture effect into account, LZ tun- 
neling seems to be unlikely, since it is exponentially sup- 
pressed for a finite 5q. When we apply the parameters of 
our experiment, the estimated LZ tunneling probability 
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Plz is negligible already for small values of Sq = 0.1: 



exp 



2hE 



exp 



AahfE,^o 



«1 (3) 



for Ejo « 45 /ieV, E'c ~ 110 /^eV and an intermediate 
amplitude a = 0.5 of rf-oscillations. On the other hand, 
experimentally we still observe a strong contribution of 
the excited band at 6q — 0.1. 

A further argument against LZ tunneling is the be- 
havior of the curves measured at ip^c = and at values 
of q around e (except for those at q — which will be 
discussed later). They clearly show the effect of qubit 
excitation, whereas the qubit energy level spacing is, in 
this case, much larger than in the case of (pdc = 

Following a model given by Ref.[2l|, we relate the 
observed mixed state to a strong coupling of non- 
equilibrium QPs to the qubit system, which allows the 
transfer of energy to the qubit and thus an excitation 
of its upper state. This QP pumping of the qubit can 
be considered as a cycle in which an unpaired non- 
equilibrium QP tunnels onto the island while the qubit 
is in the ground state |-|-), see Fig. [2] (c). As soon as 
the QP enters the transistor island, it changes the excess 
charge and induces in this way an instantaneous change 
of the working point and of the energy level splitting. 
Only when the QP tunnels to a lower-energetic state of 
the electrode, a transfer of energy to the qubit system 
occurs, exciting it to the upper state |— ), whereas a QP 
tunneling back to the initial state of the electrode does 
not induce any excitation of the qubit. According to Ref. 
[2H , the former process should prevail due to the greater 
density of the states that are close to the edge of the 
QP band in the energy spectrum [2^. Of course, such 
non-equilibrium QPs should have an energy which is at 
least by the value of Ej larger than the value of the su- 
perconductor energy gap, in order to transfer energy to 
the qubit. This would mean, however, that a QP having 
a lower energy could enter the island as well, but in that 
case the QP should leave the island without exciting the 
qubit. Trapping the QP in the island is only possible 
if the superconductor energy gap of the island is lower 
than that of the electrodes. Due to particularities of our 
sample fabrication, the energy gaps of the electrodes and 
of the transistor island are almost equal. Moreover, we 
presume that the gap of the island is even slightly larger. 
We assume that QPs having an energy lower than « Ej 
above the energy gap are available in the outer electrodes 
and that their relaxation to the gap edge occurs both via 
interaction with the lattice of the electrodes and via the 
traveling onto the island and back into the electrodes 
with simultaneous excitation of the qubit. 

The pumping process depends strongly on the qubit 
bias and on the amplitude of the rf-oscillations. This 
is demonstrated in the gate modulation curves a(VG) in 
Fig. 3 for different values of the dc phase bias. Interest- 
ingly, when increasing the amplitude of the rf-oscillations 
for ifidc — (top part of the plot) from 0.28 rad to 1.12 
rad, two dips appear around the sharp peak a,t q = e. 
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FIG. 3: Periodic gate modulation curves a{VG) for differ- 
ent values of the amplitude of the rf-oscillations, measured at 
two values of the dc phase bias. For the purpose of clarity, 
the curves are vertically shifted. The upper horizontal axis 
presents the values of the island charge. 



This is in contrast to the behavior we have described for 
Vdc = TT (bottom part of the plot). Here, the gate modu- 
lation curves practically do not depend on the amplitude 
a of the flux oscillations. We interpret these observations 
in the following way: for tp^c — 0, the pumping process is 
not active for small rf-oscillations. It starts to develop at 
increasing values of a, and the shape of the curves shows 
a strong dependence of this mechanism on the charge 
q. In particular, the pumping is partially suppressed at 
q — e even for rather large values of a. Therefore, we 
extend the theory of Ref. [2l| by establishing a selection 
rule for the QP transitions between the ground state and 
the excited state for certain flux-bias- values. 

The QP transitions are described by the tunneling 
Hamiltonian, Ht, which consists of two Hermitian con- 
jugate terms and is considered to be a perturbation term, 
see Ref. ^l|. The Josephson coupling term describ- 
ing the tunneling of Cooper pairs is naturally included 
in the Hamiltonian of the qubit. To describe this pe- 
culiar behavior observed in the experiment, we should 
take the interference effect into account which occurs 
when QPs tunnel onto the island. Let us assume that 
the energies of the non-equilibrium QP in the electrode 
Ep = (A2 + ^2)1/2 and in the island Ek = (A^ + ^2)i/2 
exceed the gap value A only slightly, i.e. the correspond- 
ing electron energies are sufficiently small, |^p|, \^k \ ^ A. 
In this case the absolute values of the Bogoliubov-Valatin 
coherence factors of the corresponding single QP states 
[13] are almost equal to each other. 



0.5. 



(4) 



This means that QPs can tunnel onto the island both 
as an electron-like particle and as a hole-like particle. 
The qubit biased by the gate charge q = {N + l)e with 
N = 0,±2,... switches into the odd (ground) state |o) 
following two paths, |±) \N) and |±) \N + 2), 



5 



as shown in Fig. [4] (a) by arrows. The sign "+" ("— ") 
stands for the even ground state (excited state) of the 
qubit. These states represent coherent superpositions of 
states with a different number of Cooper pairs on the 
island (see, e.g., Eq. (5) in Ref.[l3|)- These phase factors 
in the series do not only depend on the particular state 
of the qubit, but also on the value of the magnetic flux <i> 
which penetrates the qubit loop and thus determines the 
overall Josephson phase if — ipi + Lp2. This dependence 
is expressed most simply in the case of the small ratio 
Ej/E^, i.e., 

\±)^{\N)±e'^'^\N + 2))/^/2, (5) 

while in the general case a low- weight admixture of higher 
order terms with corresponding signs and phases, i.e., 
(±l)k^ikv/2\N + 2fc), where integer k = -1, -2, -3, . . . 
and k = 2,3,4,..., should also be added. One can 
see that in contrast to the case of a plain Cooper pair 
box which, in the degeneracy point, has a symmetric 
(antisymmetric) wave function in the ground (excited) 
state, the symmetry property of the states |±) in the 
box of the SQUID configuration alternates with the flux 
value <&. Especially, for values which correspond to an 
even number of flux quanta, $ = m$Oi™ = 0,±2,..., 
yielding e*"^/^ = 1, the ground state (excited state) of 
the qubit is described by a symmetric (antisymmetric) 
wave function, whereas for odd numbers of flux quanta, 
m = ±1,±3..., the ground state (excited state) is anti- 
symmetric (symmetric) [28j . Therefore, at any integer 
value of flux quanta in the loop, at least one of the states 
(ground state or excited state) is necessarily symmetric 
at q = e. As we see below, this property plays an impor- 
tant role for establishing the selection rule for the rate of 
the QP tunneling which includes these states. 

The operators causing transitions in the qubit are the 
charge shift operators e^*"^/^, where (j) = {ipi — </'2)/2 
is the variable conjugate to the island charge operator 
(see, e.g., Ref . [29j ) . These terms appear in the tunneling 
Hamiltonian terms as phase factors of the creation- and 
annihilation-operators associated with the island. Phys- 
ically, these phase factors ensure a coupling between the 
tunneling particle and the environment, i.e. the charge 
degree of freedom of the qubit, cf. Ref. [2^. Applying the 
Fermi Golden rule for calculating the transition rates of 
the QP tunneling onto the island, we obtain 

V^iEp,Ek) = 27r\{Ep,±\HT\o,Ek)\^S{Ep + E± ~ Ek). 

(6) 

Taking into account the simplifying assumption Eq. ([4|. 
the matrix elements in Eq. ^ can be presented as 

\{Ep,±\HT\o,Ek)\^ 

= \tpk\'[\e±f\upf\ukf + \e'^f\vp\^\vk\' 

-e±e'j.UpUkVpVk - e'^e±UpUkVpVk] 

= |ip|e±-e'±|V4, (7) 

where the matrix elements are e± = (^±\ei'l>/^\N) and 
e'^ = {±\e-"^^'^\N + 2). Here \tpk\^ = |ip is the square 



of the absolute value of the tunnel matrix element which 
is only weakly dependent on the energies Ep and E^. 
Therefore, the net rate of the QPs tunneling onto the 
island can be expressed as FiJ} oc {Sr/RT)W±, with W± — 
\e± — e'j_p (see the plots in Fig. [4] (b)) and 6r being the 
QP level spacing in the island. 

Since a.t q = e one of the qubit states |-|-) (for even 
m) or |— ) (for odd m) is described by a symmetric wave 
function, the matrix elements e+ and (for even m) 
or e_ and e'_ (for odd m) are identical, which leads to 
zero value for W+ or . This manifests itself in the zero 
minimums in the two lower curves in Fig.[4](b), as well as 
in the suppression of the corresponding transition rates 
or F'i\ respectively. In contrast to this, the transi- 
tions in which asymmetric qubit states are involved have 
an increased probability (the two upper curves in Fig.[3] 
(b)) due to the constructive interference which results 
from the opposite signs of the corresponding matrix ele- 
ments. The effect of the complete destructive interference 
occuring in transitions involving the symmetric states of 
the qubit is formulated as a selection rule. Note that the 
same selection rule is applied to transitions which are as- 
sociated with the tunneling of the QP from the island 
back into the electrode with simultaneous switching of 
the qubit into a state with symmetric wave function. If 
finite QP energies £,q and above the superconductor 
gap are taken into account, this leads to a somewhat un- 
balanced ratio of the coherence factors u and v of Eq. (4) 
and, hence, to an incomplete destructive interference by 
which the selection rule is violated. 

The destructive interference described is supposed to 
have a suppressing effect on the cyclic qubit excitation. 
In our scenario of the qubit pumping at ip^c = tt, the 
rates of both intraband and interband transitions were 
assumed to be much larger than the rate of the qubit 
relaxation B A. In this case the pumping cycle is 
closed and the two qubit states A and B are almost 
equally populated at g = e, see Fig. 2 (b). In this case, 
the average phase angle a is positioned close to, but 
slightly higher than the level of the state C, manifesting 
itself in the experimental a — Vc-curve in Fig. 2 (a) as an 
overshooting around q = e. This process is illustrated by 
the time-trace of repeated qubit transitions in the bot- 
tom panel of Fig. 4 (c) for (p^c = tt, i.e. for the case where 
neither of the states is symmetric. Let us consider now 
the phase bias cpdc — 0, where - as has been discussed be- 
fore - the qubit ground-state has a symmetric wave func- 
tion with a corresponding low rate 7symm for transitions 
into this state. Due to the effect of the destructive inter- 
ference (denoted by the dashed arrows in Fig. 1 (b)and 
Fig. 4 (a)), the rate of intraband transitions A' ^ C" is 
presumably much slower than the qubit relaxation. As a 
result, the qubit remains in the even symmetric state A' 
for a sufficiently long time or switches fast between the 
even antisymmetric state B' and the odd state C on an 
average time scale (see the time-trace for ip^c = in 
the central panel of Fig. 4 (c)). If, on the other hand, the 
flux bias is odd, e.g. (p^c = 2tt, the excited state of the 



6 



|N+2) 




(c) 

B 



).4 0.6 0.8 1.0 1.2 1.4 1.6 

q/e 



(Pdc=0 



C" 
B" 



time 



FIG. 4: (a) Energy diagram of the qubit, displaying the inter- 
fering paths of the excess QPs tunneling either - as a quasi- 
electron (to the left) - or - as a quasi-hole (to the right), (b) 
The rate factors for the intraband transition probability W+ 
and interband transition probability W~, calculated for the 
experimental parameters of the circuit for the integer number 
m of flux quanta in the loop as a function of the gate charge. 
For even m, the destructive (constructive) interference takes 
place for the intraband (interband) transitions, and for odd m 
vice versa, (c) Schematic time-traces of the cyclic transitions 
between qubit states at q — e for different values of ipdc shown 
in Fig. 1 (b). The thin arrows indicate the qubit relaxation in 
contrast to the quasiparticle-induced transitions. 



qubit is symmetric, which leads to a strong reduction 
in the interband transitions B" C". The schematic 
time-trace in the central panel in Fig. 4 (c) shows for this 
case a long series of rapid transitions A" C", inter- 
rupted by rare (with a rate of ~ 7symm) transitions into 
the state B". The escape from the excited state B" with 
an average occupation time 7"^ mainly occurs via qubit 
relaxation B" A". When we summarize the case of 



an integer dc phase bias ipdc = 27rm, the sufficiently high 
rate of relaxation with respect to the small rate due to 
the destructive interference (7^ ^ 7symm) a.t q = e en- 
sures that (a) is close to the value corresponding to the 
ground state A' or A" . Hence, this average value of the 
phase angle shows that the destructive interference effect 
suppresses the pumping cycle. 



This interference effect can qualitatively explain the 
shape of the a — Vg dependencies measured at zero flux 
value and shown in the upper part of Fig. [31 The two 
upper curves, measured at a small amplitude a of rf- 
oscillations, exhibit a shape which corresponds almost 
entirely to the expected form of the ground state. This 
behavior is related, first of all, to relatively small devi- 
ations of the Josephson Lp from the zero value, in the 
vicinity of which the destructive interference is strongest 
and, secondly, to a relatively large QP energy (~ 2E}q) 
which is required for the excitation of the qubit. The 
increase in the amplitude a leads to distinct deviations 
of the Josephson phase ip from the zero value, so the 
destructive interference is still efficient only in a narrow 
region around the value q = e. This effect manifests it- 
self in a sharp peak centered around the value q — e, 
whereas slight deviations from this value lead to the two 
side dips and are an indication of a significant admix- 
ture of the excited state. The described behavior can 
be understood as a significant weakening of the destruc- 
tive part of the interference, which is due to the devia- 
tion from the optimum point both in the dc phase and 
charge. Moreover, the excitation qubit energy is in this 
case somewhat lower than ~ 2i?jo. Note that in the case 
of a semi-integer flux quantum bias (the respective curves 
are presented in the bottom part of Fig. [3] for (^dc = 
the effect of destructive interference does not show. This 
is understandable, because neither of the wave functions 
describing the ground state and the excited state pos- 
sesses the property of symmetry (see Eq. (5) in which 
the parameter is ip = tt). 



In conclusion, we found the effect of mixing of the 
qubit states, which we explained by an energy trans- 
fer between the non-equilibrium quasiparticles and the 
two-level qubit system. To explain the qubit excitation 
observed in the experiment, we applied a picture of the 
stochastic tunneling of single quasiparticles and derived a 
selection rule for the quasiparticle tunneling at the magic 
points (7 = e of the qubit. This selection rule is applica- 
ble to transitions in which qubit states with a symmetric 
wave function are involved, and it strongly suppresses the 
mechanism of the qubit pumping. Our set-up only per- 
mitted a continuous readout of the - therefore - averaged 
qubit state. A pulsed readout set-up with a higher drive 
frequency could help to discriminate the single tunnel- 
ing events and allow probing the ground state and the 
excited state of the qubit separately. 
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